Search for: All records

A spinning black hole (BH) accreting from a disk of strongly magnetized plasma via a magnetically arrested disk is known to produce an efficient electromagnetic jet powered by the BH’s spin energy. We present general relativistic radiative magnetohydrodynamic simulations of magnetically arrested systems covering a range of sub- to super-Eddington accretion rates. Using the numerical results from these simulations, we develop formulae to describe the magnetization, jet efficiency, and spin evolution of an accreting BH as a function of its spin and accretion rate. A BH with near-Eddington accretion experiences a mild degree of spin-down because of angular momentum loss through the jet, leading to an equilibrium spin of 0.8 rather than 1.0 at the Eddington limit. As the accretion rate increases above Eddington, the spin-down effect becomes progressively stronger, ultimately converging on previous predictions based on nonradiative simulations. In particular, spin evolution drives highly super-Eddington systems toward a BH spin near zero. The formulae developed in this letter may be applied to galaxy- and cosmological-scale simulations that include BHs. If magnetically arrested disk accretion is common among supermassive BHs, the present results have broad implications for active galactic nucleus feedback and cosmological spin evolution.

 

Horizon-scale observations of the jetted active galactic nucleus M87 are compared with simulations spanning a broad range of dissipation mechanisms and plasma content in three-dimensional general relativistic flows around spinning black holes. Observations of synchrotron radiation from radio to X-ray frequencies can be compared with simulations by adding prescriptions specifying the relativistic electron-plus-positron distribution function and associated radiative transfer coefficients. A suite of time-varying simulations with various spins, plasma magnetizations and turbulent heating and equipartition-based emission prescriptions (and piecewise combinations thereof) is chosen to represent distinct possibilities for the M87 jet/accretion flow/black hole system. Simulation jet morphology, polarization, and variation are then ‘observed’ and compared with real observations to infer the rules that govern the polarized emissivity. Our models support several possible spin/emission model/plasma composition combinations supplying the jet in M87, whose black hole shadow has been observed down to the photon ring at 230 GHz by the Event Horizon Telescope (EHT). Net linear polarization and circular polarization constraints favour magnetically arrested disc (MAD) models whereas resolved linear polarization favours standard and normal evolution (SANE) in our parameter space. We also show that some MAD cases dominated by intrinsic circular polarization have near-linear V/I dependence on un-paired electron or positron content while SANE polarization exhibits markedly greater positron-dependent Faraday effects – future probes of the SANE/MAD dichotomy and plasma content with the EHT. This is the second work in a series also applying the ‘observing’ simulations methodology to near-horizon regions of supermassive black holes in Sgr A* and 3C 279.

 

We present two general relativistic radiation magnetohydrodynamics (GRRMHD) simulations of magnetically arrested discs (MADs) around non-spinning (a* = 0) and spinning (a* = 0.9) supermassive black holes (BHs). In each simulation, the mass accretion rate is decreased with time such that we sample Eddington-scaled rates over the range $3 \gtrsim \dot{M}/\dot{M}_{\rm {Edd}}\gtrsim 0.3$. For the non-spinning BH model, the total and radiative efficiencies increase as the accretion rate decreases, varying over the range $\eta _{\rm {tot}}\sim 9\!-\!16{{\ \rm per\ cent}}$ and $\eta _{\rm {rad}}\sim 6{-}12{{\ \rm per\ cent}}$, respectively. This model shows very little jet activity. In contrast, the spinning BH model has a strong relativistic jet powered by spin energy extracted from the BH. The jet power declines with accretion rate such that $\eta _{\rm {jet}}\sim 18{-}39{{\ \rm per\ cent}}$ while the total and radiative efficiencies are $\eta _{\rm {tot}}\sim 64{-}100{{\ \rm per\ cent}}$ and $\eta _{\rm {rad}}\sim 45{-}79{{\ \rm per\ cent}}$, respectively. We confirm that mildly sub-Eddington discs can extract substantial power from a spinning BH, provided they are in the MAD state. The jet profile out to $100\, GM/c^2$ is roughly parabolic with a power-law index of k ≈ 0.43−0.53 during the sub-Eddington evolution. Both models show significant variability in the outgoing radiation which is likely associated with episodes of magnetic flux eruptions. The a* = 0.9 model shows semiregular variations with a period of $\sim 2000\, GM/c^3$ over the final $\sim 10\, 000\, GM/c^3$ of the simulation, which suggests that magnetic flux eruptions may be an important source of quasi-periodic variability. For the simulated accretion rates, the a* = 0 model is spinning up while the a* = 0.9 model is spinning down. Spinup–spindown equilibrium of the BH will likely be achieved at 0.5 < a*, eq < 0.6, assuming continuous accretion in the MAD state.

 

Tidal disruption events (TDEs) around supermassive black holes (SMBHs) are a potential laboratory to study super-Eddington accretion disks and sometimes result in powerful jets or outflows which may shine in the radio and sub-millimeter bands. In this work, we modeled the thermal synchrotron emission of jets by general relativistic radiation magneto-hydrodynamics (GRRMHD) simulations of a BH accretion disk/jet system which assumed the TDE resulted in a magnetized accretion disk around a BH accreting at ∼12–25 times the Eddington accretion rate. Through synthetic observations with the Next Generation Event Horizon Telescope (ngEHT) and an image reconstruction analysis, we demonstrate that TDE jets may provide compelling targets within the context of the models explored in this work. In particular, we found that jets launched by a SANE super-Eddington disk around a spin a*=0.9 reach the ngEHT detection threshold at large distances (up to 100 Mpc in this work). A two-temperature plasma in the jet or weaker jets, such as a spin a*=0 model, requires a much closer distance, as we demonstrate detection at 10 Mpc for limiting cases of a*=0,R=1 or a*=0.9,R=20. We also demonstrate that TDE jets may appear as superluminal sources if the BH is rapidly rotating and the jet is viewed nearly face on. 

We present general relativistic radiation magnetohydrodynamics (GRRMHD) simulations of super-Eddington accretion flows around supermassive black holes (SMBHs), which may apply to tidal disruption events (TDEs). We perform long duration ($t\ge 81,200\, GM/c^3$) simulations that achieve mass accretion rates ≳11 times the Eddington rate and produce thermal synchrotron spectra and images of their jets. Gas flowing beyond the funnel wall expands conically and drives a strong shock at the jet head while variable mass ejection and recollimation, along the jet axis, results in internal shocks and dissipation. Assuming the ion temperature (Ti) and electron temperature (Te) in the plasma are identical, the radio/submillimetre spectra peak at >100 GHz and the luminosity increases with BH spin, exceeding $\sim 10^{41} \, \rm {erg\, s^{-1}}$ in the brightest models. The emission is extremely sensitive to Ti/Te as some models show an order-of-magnitude decrease in the peak frequency and up to four orders-of-magnitude decline in their radio/submillimetre luminosity as Ti/Te approaches 20. Assuming a maximum VLBI baseline distance of 10 Gλ, 230 GHz images of Ti/Te = 1 models shows that the jet head may be bright enough for its motion to be captured with the EHT (ngEHT) at D ≲ 110 (180) Mpc at the 5σ significance level. Resolving emission from internal shocks requires D ≲ 45 Mpc for both the EHT or ngEHT.

 

We explore the plasma matter content in the innermost accretion disk/jet in M87* as relevant for an enthusiastic search for the signatures of anti-matter in the next generation of the Event Horizon Telescope (ngEHT). We model the impact of non-zero positron-to-electron ratio using different emission models, including a constant electron to magnetic pressure (constant βe model) with a population of non-thermal electrons as well as an R-beta model populated with thermal electrons. In the former case, we pick a semi-analytic fit to the force-free region of a general relativistic magnetohydrodynamic (GRMHD) simulation, while in the latter case, we analyze the GRMHD simulations directly. In both cases, positrons are being added at the post-processing level. We generate polarized images and spectra for some of these models and find out that at the radio frequencies, both of the linear and the circular polarizations are enhanced with every pair added. On the contrary, we show that, at higher frequencies, a substantial positron fraction washes out the circular polarization. We report strong degeneracies between different emission models and the positron fraction, though our non-thermal models show more sensitivities to the pair fraction than the thermal models. We conclude that a large theoretical image library is indeed required to fully understand the trends probed in this study, and to place them in the context of a large set of parameters which also affect polarimetric images, such as magnetic field strength, black hole spin, and detailed aspects of the electron temperature and the distribution function. 

ABSTRACT We use the general relativistic radiation magnetohydrodynamics code KORAL to simulate the accretion disc formation resulting from the tidal disruption of a solar mass star around a supermassive black hole (BH) of mass 106 M⊙. We simulate the disruption of artificially more bound stars with orbital eccentricity e ≤ 0.99 (compared to the more realistic case of parabolic orbits with e = 1) on close orbits with impact parameter β ≥ 3. We use a novel method of injecting the tidal stream into the domain, and we begin the stream injection at the peak fallback rate in this study. For two simulations, we choose e = 0.99 and inject mass at a rate that is similar to parabolic TDEs. We find that the disc only becomes mildly circularized with eccentricity e ≈ 0.6 within the 3.5 d that we simulate. The rate of circularization is faster for pericenter radii that come closer to the BH. The emitted radiation is mildly super-Eddington with $L_{\rm {bol}}\approx 3{-}5\, L_{\rm {Edd}}$ and the photosphere is highly asymmetric with the photosphere being significantly closer to the inner accretion disc for viewing angles near pericenter. We find that soft X-ray radiation with Trad ≈ 3–5 × 105 K may be visible for chance viewing angles. Our simulations suggest that TDEs should be radiatively inefficient with η ≈ 0.009–0.014. 

We present the results of nine simulations of radiatively inefficient magnetically arrested discs (MADs) across different values of the black hole spin parameter a*: −0.9, −0.7, −0.5, −0.3, 0, 0.3, 0.5, 0.7, and 0.9. Each simulation was run up to $t \gtrsim 100\, 000\, GM/c^3$ to ensure disc inflow equilibrium out to large radii. We find that the saturated magnetic flux level, and consequently also jet power, of MAD discs depends strongly on the black hole spin, confirming previous results. Prograde discs saturate at a much higher relative magnetic flux and have more powerful jets than their retrograde counterparts. MADs with spinning black holes naturally launch jets with generalized parabolic profiles whose widths vary as a power of distance from the black hole. For distances up to 100GM/c2, the power-law index is k ≈ 0.27–0.42. There is a strong correlation between the disc–jet geometry and the dimensionless magnetic flux, resulting in prograde systems displaying thinner equatorial accretion flows near the black hole and wider jets, compared to retrograde systems. Prograde and retrograde MADs also exhibit different trends in disc variability: accretion rate variability increases with increasing spin for a* > 0 and remains almost constant for a* ≲ 0, while magnetic flux variability shows the opposite trend. Jets in the MAD state remove more angular momentum from black holes than is accreted, effectively spinning down the black hole. If powerful jets from MAD systems in Nature are persistent, this loss of angular momentum will notably reduce the black hole spin over cosmic time.